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ABSTRACT 


In  this  paper,  we  present  an  effective  technique  to  detect  and  localize  a moving  target  behind  the  wall.  The 
method  based  on  one  transmitter  and  two  receiver  antenna  system  configuration.  The  precision  of  the  method  is  first 
surveyed  by  simulating  a hypothetical  target.  Then,  the  algorithm  is  tested  with  a real  data  gathered  from  a measurement 
that  was  accomplished  behind  a brick  wall  with  an  ultra- wide-band  (UWB)  transceiver  system.  Thru-the-wall  radar 


imaging  results  show  that  algorithm  successfully  detects  and  indicates  the  correct  locations  of  a moving  person  behind  the 
wall.  Our  numerical  and  experimental  studies  in  this  work  clearly  demonstrate  that  the  proposed  method  can  be 


appropriately  used  in  similar  TWR  applications  or  localization  problems. 

KEYWORDS:  Remote  Sensing,  Through-The-Wall  Radar,  Moving  Target  Detection,  Radar  Signal  Processing, 
Localization  Algorithm 

1.  INTRODUCTION 

Through-wall  radar  (TWR)  is  has  been  an  highly  developing  technology  for  detecting  and  localizing 
persons/objects  behind  the  wall  or  any  opaque  obstacle  [1-3].  TWR  has  many  practices  ranging  from  civil  applications 
such  as  earthquake/fire  rescue  operations  to  military  ones  including  detecting  terrorists  inside  buildings  [1-5].  For  the  TWR 
design,  ultra- wide  band  (UWB)  time-domain  radar  systems  are  thought  to  be  more  advantageous  to  frequency  domain 
counterparts  since  they  can  operate  faster  and  also  maintain  high  resolution  features  with  much  higher  signal  power  levels 
for  better  penetration  and  deeper  detection  ranges  [5-6].  For  these  advantageous  properties,  many  scientists  have  utilized 
UWB  time-domain  TWR  usage  by  applying  various  detection/localization  techniques  to  focus  and  image  objects  or  human 
movement  behind  walls  [6-10].  One  of  the  promising  methods  is  called  the  back-projection  algorithm  (BPA)  which  has 
been  also  applied  successfully  in  conventional  synthetic  aperture  radar  (SAR)  applications  [11-14].  Although  BPA  has  been 
frequently  used  in  TWR  imaging  usages,  its  summation  step  causes  undesired  artifacts  in  the  form  of  long  side  lobes  that 
negatively  affect  the  resolutions  in  the  constructed  images.  Similar  to  BPA,  other  focusing  techniques  such  as  Kirchoff 
Migration  [15]  and  f-k  migration  [16]  have  the  main  disadvantage  of  being  slow  due  to  iterative  summation  step  within  the 
algorithm.  In  addition  to  these  well-known  robust  focusing  methods,  some  high-resolution  techniques  have  also  been 
applied  by  some  researchers  in  pursuit  of  having  better  detection  outcomes  [17-19].  One  of  the  most  promising  candidate  is 
the  beams  pace  multiple  signal  classification  algorithm  (BS -MUSIC)  that  have  the  ability  to  provide  highly  resolved 
images  [17].  On  the  other  hand,  it  requires  the  a priori  knowledge  of  some  parameters  of  scene  such  as  wall  thickness  and 
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the  wall’s  dielectric  constant.  Furthermore,  MUSIC  is  quite  slow  while  computing  the  estimate  of  the  covariance  matrix 
and  array  response  vector.  Sparse  Tomographic  Inverse  Scattering  Approach  [18]  and  an  autofocusing  technique,  based  on 
higher  order  statisticsmethod  [19]  have  also  been  applied  for  seeking  a good  performance  on  detection  with  a fine 
resolution  metrics.  They  require  coherent  processing  of  the  collected  data  sets  that  cannot  be  practical  in  most  TWR 
applications. 

In  this  work,  we  have  presented  a non-coherent  technique  based  on  geometrical  bilateration  approach  for 
one-transmitter  two-receiver  TWR  set-up.  This  method  mainly  calculates  the  intersection  point  of  the  two  ellipses  that  are 
formed  by  one-transmitter  one-receiver  pairs.  This  paper  is  organized  as  follows:  In  the  next  section,  the  imaging  algorithm 
based  on  intersection  of  ellipses  is  described  with  details.  In  Section  3,  the  simulation  results  are  presented  to  assess  the 
performance  of  the  algorithm.  In  Section  4,  the  algorithm  is  tested  by  the  real  data  collected  with  a UWB  time  domain 
radar  system  to  detect  and  locate  a moving  human  behind  a brick  wall.  In  the  last  section,  the  discussion  about  the  cons  and 
pros  of  the  implemented  algorithm  are  given  and  the  work  is  concluded. 

2.  LOCALIZATION  ALGORITHMS ASED  ON  BILATERATION  TECHNIQUE 

Our  detection  and  localization  algorithm  relies  on  the  assumptions  that  the  wall  thickness  is  small  compared  to 
target  range  and  the  dielectric  permittivity  of  the  wall  medium  ( er ) is  not  high.  In  fact,  these  assumptions  generally  hold 
true  provided  that  the  practical  wall  thickness  generally  ranges  from  15  cm  to  40  cm  and  er  varies  from  1.5  to  6 for 
wooden,  brick  and/or  concrete  walls  [20].  With  this  construct,  the  amount  of  shift  between  the  actual  location  and  the 
estimated  location  of  the  target  ( dshift ) is  negligible  as  illustrated  in  Figure  1.  The  error  associated  with  this  assumption 
has  been  studied  in  detail  in  [21]. 


The  principle  of  our  algorithm  is  based  on  multi-static  geometry  in  which  a transmitting  antenna  (Tx)  and  two 
receiving  antennas  (Rx1and  Rx2)are  positioned  behind  the  wall  as  illustrated  in  Figure  1.  As  shown  from  the  figure,  as  the 
Tx  is  situated  in  the  middle  and  Rx^nd  Rx2  are  put  at  the  same  distance  of  D away  from  Tx  along  the  cross-range  axis  (x). 
The  distance  from  Tx  to  point  target;  r0  and  the  distances  from  point  target  to Rxt(i  = 1, 2);  rx  andr2,  are  estimated  by  the 
help  of  round  trip  time  and  the  speed  of  the  electromagnetic  (EM)  wave.  In  fact;  the  locations  of  Txand  Rxt  scorrespond  to 
the  focal  points  of  the  ellipse  curves.  Therefore;  the  sum  of  the  distances  from  the  foci  to  any  points  on  the  ellipses  (r0  + 
rlfr0  + r2)  are  constant  and  they  are  equal  to  the  horizontal  major  axis  of  2 a.  For  our  radar  geometry;  two  elliptic  curves 
are  occurred  due  to  two  receivers.  The  equations  of  these  two  ellipses  centered  at  (—D/2,0)  and  (D/2,0)  with  respect  to 
the  x and  y-axis  as; 


C x+D/2 )2  y2  _ 

ax2  + b2 

(x-U/2)2  y2  _ 
a22  "h  b2 


(1) 

(2) 


Here;  semi-major  axes  (a1anda2);  and  semi-minor  axes  (h1andh2)  of  two  ellipses  according  to  Figure  1 are 
described  as  below: 


ai  = r0  + r± 
a2  =r0+  r2 
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^1  — V ai2  — (D/2)2 

b2  = V«22  - (£»/2)2 


It  is  assumed  that  there  is  a perfect  point  scatter  as  a target  behind  the  wall  as  shown  in  Figure.  1 and  its  location 
F(x0,y0)  can  be  simply  obtained  by  intersection  of  these  two  ellipses  which  are  defined  in  Eqn.  (1)  and  Eqn.  (2), 
respectively.  The  formulation  steps  in  calculating  the  intersection  of  these  curves  are  expressed  briefly  as  follows; 

A = a22  ■ bt2  - at2  ■ b22 

B = D ■ (a22  ■ bx2  + ■ b2)  (5) 

C = (D/2)2  ■ A + at2  ■ a22{b2 2 - bt2) 


*1,2 


B+y/B2-4AC 

2A 


y = 


(6) 


Wherei4,  B and  C are  the  coefficients  of  the  quadratic  polynomial  of  Ax2  + Bx  + C = 0.  Then,  the  cross-range 
location  (x)  and  the  range  position  (y)  of  the  point  target  behind  the  wall  can  be  calculated  by  finding  the  roots  of  this 
polynomial  and  solving  Eqn.  (6).  In  fact,  this  equation  set  yields  two  dual  solutions  of  ( x1,y1 ) and  (x2,y2).  Based  on  the 
assumed  geometry  in  Figure  1 ; of  course,  the  positive  notation  for  the  range  value  (y)  should  be  chosen  as  the  solution  that 
will  give  only  one  set  of  either  (xlf  yx)  or  (x2,  y2). 


SIMULATION  SCENARIO 


We  have  constructed  a simulation  scenario  in  Matlab  programming  [22]  to  assess  the  performance  of  the 
proposed  detection  and  localization  algorithm.  For  this  goal,  a hypothetical  moving  target  was  considered.  In  the 
simulation,  we  have  chosen  the  separation  of  antennas  atD  = 30  cmas  seen  in  Figure  2 where  the  geometry  of  the  scenario 
is  given.  The  frequency  of  the  radar  system  was  altered  from  1 GHz  to  6GHz  for  a total  of  200  discrete  frequencies  in 
stepped-frequency  radar  configuration.  Firstly,  the  starting  point  of  moving  target  was  assumed  at  point/5 (1,3)  on  the 
range  (y)and  cross-range  (x)plane.  Then,  it  was  reached  up  to  point  P ^_1,3^bya  curved  path  for  a total  of  18  distinct  target 
locations  as  shown  in  Figure  2.  So,  the  frequency  diverse  two-dimensional  (2D)  scattered  electric  field  data;  Ff(/c,  r) and 
Ff  (/c,  r)were  collected  by  the  receivers  at  different  time  instants  of  r^as 


Ei(k,Ti)  = CTj  ■ e Mroi+r-it)  . = ^ ^ ^ 

E2(k,Ti)  = a2  ■ e-,k(ro;+n>i) 


(10) 


Here;  a1  and  <j2  are  the  EM  reflectivity  constants  for  the  first  and  the  second  receiver,  respectively.  For  this 
simulation,  they  were  assumed  to  be  independent  from  frequency  and  aspect  angle  so  that  ox  = cr2  = 1 . Now,  we  can  get 
the  bi-static  ranges  of  these  collected  scattered  electric  field  data  by  taking  the  one  dimensional  (ID)  inverse  Fourier 
transform  (IFT)  of  Ff(/e,  r)  and  Ff(/c,  r)  with  respect  to  wave-number;  k.  Then,  the  2D  range-time  scattered  field  vector 
sets  of  E^(r,z)  and  F|(r,  r)  were  obtained  as  plotted  in  Figure  3 (a)  and  Figure  3 (b),  respectively.  As  shown  from  these 
range  profiles,  the  movement  of  point  scatterer  can  be  easily  observed.  Next,  the  terms  of  semi-major  axes (ax and  a2)and 
semi-minor  axes  and  b2)  were  calculated  that  were  explained  in  the  previous  section.  Then,  the  steps  of  the  algorithm 
were  employed  in  the  above-mentioned  successive  order  to  obtain  the  range  (y)  and  cross-range  (x) coordinates  of  18 
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different  scene  snapshots.  In  Figure  4,  employment  of  the  bilateration  technique  has  been  realized  for  three  different  time 
instants.  The  ellipse  that  corresponds  to  the  first  receiver  was  drawn  as  red  dashed  line  whereas  the  one  that  corresponds  to 
the  second  receiver  was  drawn  as  blue  dashed  line.  The  intersections  of  the  ellipses  are  represented  as  black  dots.  In 
Figure  5,  both  the  real  locations  (drawn  as  black  circles)  and  the  estimated  locations  (plotted  as  magenta  crosses)  by  our 
algorithm  of  the  moving  target  for  different  1 8 time  instants  are  plotted.  The  visual  comparison  of  these  two  results  clearly 
demonstrates  how  our  algorithm  successfully  detects  and  pinpoints  the  locations  of  the  moving  point-scatterer  with  a 
reliable  prediction. 

MEASUREMENT  RESULTS 

After  validating  the  accuracy  of  the  algorithm  by  the  help  of  synthetic  scenario,  the  algorithm  was  tested  with  a 
real  data  collected  by  an  experimental  test-bed.  For  this  goal,  an  experimental  scenario  whose  geometry  as  shown  in 
Figure  6 was  formed.  Tx  antenna  was  placed  at  (0,  0)  and  two  Rx  antennas  were  positioned  equally  spaced  as  30  cm  on 
either  side  of  it.  Additionally;  the  antenna  configuration  was  positioned  5 cm  away  against  a 20  cm  thick  brick  wall. 
During  the  experiment,  a male  human  was  stood  behind  the  wall  and  walked  away  and  towards  the  wall  twice  during  the 
time  period  of  21  5.  As  it  can  be  observed  clearly  from  the  Figure  6;  the  man  started  to  walk  at  a radial  distance  of  0.35  m 
from  the  Tx  antenna  and  returned  from4.1  m in  the  first-round  trip.  In  his  second-round  walk;  he  followed  the  same  route 
and  trip  was  ended  at  2.4  m radial  distance  from  the  transmitter. 

Our  experimental  measurement  system  is  shown  in  Figure  7.  The  measurement  set-up  was  composed  of  an 
Analog/Digital  convertor  (Geozondas  GZ6E)  and  pulse  generator  head  (GZ1118GN-01EV)  which  was  used  as  the  time 
domain  transceiver  system  as  seen  Figure  7 (a).  Also,  we  used  one  Tx  and  two  Rx  horn  antennas  that  can  work  between  1 
and  6 GHz.  For  the  automatic  control  of  the  whole  system,  we  have  developed  a Matlab  script.  While  the  man  was  walking 
behind  the  wall  (see  Figure  7 (b)),  we  collected  consecutive92time  domain  responses  of  the  scene  during  the  21  5 period. 
The  collected  raw  data  from  Rx i and  Rx2  are  as  depicted  in  Figure  8 (a)  and  (b),  respectively.  We  can  see  from  the  figures 
that  signals  which  were  shown  at  r = 0 represented  the  cross-talk  between  Tx  and  the  Rx  antennas  and  also  taken  as  the 
phase  center  of  our  experiment.  The  wall  was  seen  from  r = 0.15  m to  r = 0.35  m as  the  highest  echo  shown  up  as  the 
straight  lines  in  all  ID  range  profiles.  Then,  we  have  applied  the  background  subtraction  procedure  to  raw  dataset  to 
improve  the  signal-to-noise  ratio(SNR)  and  also  to  remove  the  artifacts  caused  by  the  wall.  Sequential  differentiation 
process  has  been  applied  by  subtracting  the  images  in  a sequential  order.  Thus,  the  stationary  objects  and  other  undesirable 
clutters  were  suppressed  and  the  reflections  of  man’s  walking  were  substantially  raised  up  as  shown  in  Figure  9 (a)  and  (b) 
for  the  first  and  the  second  receiver,  respectively.  When  t = 5.8  s,  the  man  was  about  3.75  m away  from  the  wall  and  it 
took  6.1  5 to  reach  the  nearest  position  and  he  returned  the  same  range  location  as  shown  in  Figure  9.  Weapplied  the  steps 
of  our  implementation  for  the  bilateration  technique  that  was  described  in  Section  2.  Figure  10  shows  the  red  and  blue 
elliptic  curves  that  were  collected  by  Rx i and  Rx2 , respectively.  The  black  dots  on  these  curves  indicate  the  intersection 
point  of  these  two  ellipses  calculated  with  the  proposed  method  at  different  time  instants.  In  Figure  11,  the  estimated 
locations  of  the  person  between  t = 5.8  s and  t = 17.1  s are  plotted.  These  results  are  in  very  good  agreement  with  the  real 
positions  of  the  walking  target  during  his  round-trip  movement. 

CONCLUSIONS 

In  this  paper,  we  have  proposed  and  implemented  a method  based  on  the  range  profile  knowledge  from  the 
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receiving  antennas  that  was  positioned  at  both  side  of  the  transmitting  antenna  for  TWR  usage  to  detect  and  localize  a 
target  behind  an  opaque  obstacle.  The  algorithm  relies  on  the  estimation  of  intersecting  the  ellipses  formed  by  1-Tx  and 
1-Rx  pair  for  the  two  receivers.  The  proposed  method  was  examined  with  both  a simulation  data  constructed  in  Matlab 
environment  and  the  real  data  from  a UWB  time-domain  radar  system.  The  resultant  detection  and  localization 
performance  of  the  algorithm  clearly  demonstrate  the  validity  and  the  success  of  the  proposed  approach. 

This  algorithm  has  the  advantage  of  being  very  fast  since  the  location  of  the  target  is  directly  calculated  by  taking 
ID  range  profile  measurement  and  applying  the  analytical  formula  given  in  Eqns.  (5)  and  (6).  Therefore,  real-time 
detection  and  localization  of  the  target  is  realized.  On  the  other  hand,  this  approach  has  the  disadvantage  of  detecting  only 
one  movement  since  the  intersection  line  of  more  than  two  ellipses  will  produce  many  intersections.  While  some  of  these 
intersections  (i.e.,  1-1,  2-2,  3-3,..  etc.)  would  represent  the  exact  locations  of  the  target,  other  intersecting  locations 
(i.e.,  1-2,  2-1,  1-3,  3-1  2-3,  3-2,..  etc.),  would  produce  hypothetical,  ghost  targets.  As  the  future  research;  therefore,  this 
algorithm  can  be  further  improved  to  detect  and  locate  multiple  targets  behind  the  wall  by  taking  care  of  this  ghost  crossing 
problem. 
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FIGURE  CAPTIONS 
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Figure  1:  Through-The-Wall  Radar  Geometry  of  Bilateration  Method 
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Figure  2:  Geometry  of  Simulation  Scenario 
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Figure  3:  Bistatic  Range  Profiles  for  (A)  Receiver  #1  and  (B)  Receiver  #2  for  Different  Time  Instants  from  ^To  t18 
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Figure  4:  Elliptic  Curves  Related  to  Receiver  #1  [Drawn  As  Red]  and  Receiver  #2  [Drawn  As  Blue]  and  Intersection 
Points  of  These  Curves  for  Three  Different  Time  Instants  for  Simulation  Results 
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Figure  5:  True  Locations  (Black  Circles)  and  Estimated  Locations  (Magenta  Crosses)  of  Human 

Target  Locations  after  the  Algorithm 
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Figure  6:  The  Geometry  of  the  Measurement  Geometry:  Human  Walking  Towards  and  Away  from  the  Wall 
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Figure  7:  (a)  A Picture  of  the  Experimental  Set-Up,  (b)  A Scene  From  the  Measurement:  Human  Walking 

Away  and  Towards  the  Wall 
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Figure  8:  The  Raw  Range-Time  Images  of  the  Collected  Data  Sets 
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Figure  9:  The  Differentiated  Range-Time  Images  of  the  Collected  Data  Sets 
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Figure  10:  Measured  Elliptic  Curves  Related  to  Receiver  #l[Drawn  as  Red]  and  Receiver  #2  [Drawn  as  Blue]  and 
Intersection  Points  of  Thesecurves  for  Three  Different  Time  Instants 


E 


o 

o 


o 


o 

Cj 


JO  = 5.8  s 


j t = 9.55  s 


Range  [m] 


Range  [m] 


Figure  11:  Calculated  Locations  of  the  Walking  Man  after  Applying  the  Algorithm 
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